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C . Reports  published  and  accepted  for  publication. 

The  research  projects  supported  entirely  or  in  part  by 
this  grant  are  described  in  forty-seven  research  papers  and 
books.  These  are  summarized  below  by  area  of  application. 

A.  Electromagnetic  phenomena. 

1 . Electromagnetic  wave  propagation. 

ARO-10,  "Perturbation  Theory  of  Nonlinear  Electromagnetic 
Wave  Propagation"  by  Joseph  B.  Keller  and  Martin  H.  Millman, 

The  Physical  Review,  181,  5,  1730-1747  (1969). 

A perturbation  method,  previously  employed  to  treat  nonlinear  lioundary-vatue  problems  involving 
partial  difTerential  equations,  is  used  to  study  nonlinear  electromagnetic  n-ave  propagation.  The  problems 
considered  are  I.  electromagnetic  wave  propagation  in  a waveguide  containing  a nonlinear  i.^tropic  me- 
dium; II.  electromagnetic  osdllations  in  a cavity  containing  a nonlinear  dispersive  anisotropic  meilium; 

III.  propagation  and  interaction  of  electromagnetic  waves  in  a nonlinear  dispersive  anisotropic  ntedium; 
rV.  reflection  of  an  electromagnetic  wave  from  a nonlinear  dispersive  unia.tial  medium.  The  method  avoids 
secular  terms,  which  arise  in  some  perturbation  treatments.  The  results  sliow  how  the  wave  number,  propa- 
gation velocity,  and  angle  of  refraction  of  a wave  depend  upon  its  amplitude. 

ARO-47,  "Dynamics  of  the  Josephson  Junction"  by 
Frank  C.  Hoppensteadt , Willard  L.  Miranker  and  Mark  Levi,  Quart. 
Jour,  on  Appl.  Math.,  in  press. 


2 . Nonlinear  circuit  theory. 


ARO-35,  "Asymptotic  Stability  in  Networks  with  Parasites" 
by  Frank  Hoppensteadt , Proceed.  11th  Annual  Allerton  Conf.  on 
Circuit  and  System  Theory,  703-707  (1973) . 

An  Iraportant  question  is  how  the  presence  of  parasitics  can  influence 
the  large-time  stability  of  an  electrical  network.  Many  investi [Rations  of 
this  problem  have  been  based  on  the  construction  of  Liapunov  functions,  a 
method  which  can  become  quite  tedious.  However,  frequently  much  more  in- 
formation can  be  obtained  under  the  same  conditions  by  using  a well  knovm 
perturbation  method  to  construct  approximate  solutions  of  the  problem. 
These  lead  to  useful  approximations  to  the  domain  of  stability,  to  the 
large-time  state  of  the  system,  and  in  asymptotically  steady  problems,  to 
the  system's  steady  state. 


ARO-36,  "Numerical  Solution  of  Differential  Equations  with 
Rapidly  Changing  Solutions"  by  Frank  Hoppensteadt  and 
Willard  L.  Miranker,  IBM  Tech.  Report  RC4792,  1-39  (1974). 

Initial  value  problems  for  weakly  nonlinear  systems  of  ordinary 
differential  equations  are  studied  here.  First  a perturbation  result  is 
obtained  which  is  shown  to  give  a valid  approximation  to  the  solution. 

The  leading  term  in  the  approximation  expansion  is  determined  by  solving 
two  associated  problems.  This  approximation  is  particularly  relevant  to 
problems  which  have  rapidly  oscillating  and  rapidly  decaying  solutions. 
Next,  a numerical  procedure,  based  on  the  asymptotic  approximation,  is 
formulated.  The  procedure  bypasses  a major  technical  problem  in  the 
construction  of  the  approximation,  and  it  provides  an  attractive  method 
for  computing  the  solution. 


ARO-45,  "Frequency  Entrainment  of  a Forced  Van  der  Pol 
Oscillator"  by  J.  E.  Flaherty  and  F.  C.  Hoppensteadt,  Studies 
in  Appl.  Math.,  in  press. 

A van  der  Pol  relaxation  oscillator  that  is  subjected  to  external 
sinusoidal  forcing  can  exhibit  stable  and  unstable  periodic  and  almost 
periodic  responses.  For  some  forcing  amplitudes  it  even  happens  that 
two  stable  subharmonic  having  different  periods  may  coexist.  We  in- 
vestigate here  the  stctble  responses  of  such  forced  oscillators.  By 
numerically  computing  the  rotation  number  of  stable  oscillations  for 
various  values  of  the  forcing  eunplltude  and  oscillator  tuning,  we 
obtain  descriptions  of  regions  of  phase  locking,  successive  bifurcation 
of  stable  subharmonic  and  almost  periodic  oscillations,  and  overlap 
regions  where  two  distinct  stable  oscillations  can  coexist. 


3. 


Propagation  in  random  media. 


ARO-18,  "Stochastic  Differential  Equations  with  Applica 
tions  to  Random  Harmonic  Oscillators  and  Wave  Propagation  in 
Random  Media"  by  George  Papanicolaou  and  Joseph  B.  Keller, 
SIAM  J.  Appl.  Math.,  21,  2,  287-305  (1971). 


Th^'  two-iimc  mclhoil  used  to  obtain  an  c\p.insion.  valid  for  i.  small  and  t larpi.  of 
the  vfclor  soiimon  of  an  abstract  ordinary  differential  cnualion  involviiis;  t.  Tl;c  s.mic  method  is 
used  lo  get  evpansions  of  functions  of  n.  The  results  are  shown  to  apply  to  the  solutions  of  stocliaslic 
equations.  I hey  are  used  to  find  the  first  two  moments  and  the  transition  probability  of  the  displace- 
rieiii  of  a harmonic  o.sei!laior  with  sprini:  constant  a random  function  of  r.  The  rcMiilt  com  tins  the 
condition  for  mean  square  stability  due  to  Stratonovich.  The  results  are  also  applied  to  ..iie-siiiiiensional 
wave  propagation  through  a layer  with  refractive  index  a random  function  of  position.  Tin  v are  used 
to  find  the  mean  square  amplitude  relleclioii  and  transmission  coeflicienis.  which  arc  ius!  the  mean 
power  reHeetion  and  transmission  coellicients.  A graph  of  the  mean  square  transmission  eoeilieieni  as 
a function  of  layer  thickness  is  preseniej.  The  results  arc  ul>o  compared  with  those  oSiain.ible  by  other 
methods. 


ARO-13,  "Finite-Amplitude  Sound-Wave  Propagation  in  a 


Waveguide"  Joseph  B.  Keller  and  Martin  H.  Millman,  J.  Acoust. 
Soc.  Am.,  49,  1 (Part  2),  329-333  (1971). 

A sound  wave  with  angular  frequency  u and  finite  amplitude  «,  propagating  in  a rigid  waveguide  of  any  non- 
rectangular  cross  section,  is  considered.  The  pressure  and  velocity  potential  of  the  wave,  corresponding  to 
each  mode  of  the  linear  theory,  are  determined  as  power  series  in  t.  The  leading  term  is  just  the  mode  given 
by  the  linear  theor>-  and  further  terms  are  finite-amplitude  corrections  to  it.  The  wavenumber  or  propagation 
constant  k,  of  the  nth  mode  is  found  to  depend  upon  < and  is  also  determined  as  a power  series.  The  fact  that 
depends  upon  ,,  and  the  e.xpression  for  4,(«),  is  the  most  interesting  consequence  of  the  analysis.  This 
shows  that  the  propagation  velocity  «/*,  depends  upon  e.  The  results  are  S|iecialued  to  a waveguide  of 
circular  cross  section.  Then  numerical  values  of  *.(,)  are  given  for  a circular  guide  filled  with  air  for  n- 1,  2, 

3,  4,  and  5.  The  lowest  mode  (»i=0)  is  not  considered  because  shocks  occur  in  it.  The  method  of  analysis 
is  a perturbation  e.\pansion  adapted  to  eliminate  secular  terms.  It  has  been  used  before  to  treat  periodic 
finite-amplitude  sound  waves  in  a closed  conuiner  and  periodic  finite-amplitude  vibrations  of  strings  and 
beams  by  J.  B.  Keller  and  L.  Ting  [Comm.  Pure  .\ppl.  Math.  19, 37i  (1966)],  and  to  ueat  nonlinear  electro-  ' 
magnetic  wave  propagation  [Phys.  Rev.  181, 1730  (1969)]  and  other  nonlinear  boundary-value  problems  by 
the  authors  [J.  Math.  Phys.  10,  342  (1969)].  j 


ARO-32,  "Nonlinear  Forced  and  Free  Vibrations  in  Acoustic 
Waveguides"  by  Joseph  B.  Keller,  J.  Acoust.  Soc.  Am.,  3, 
524-527  (1974) . 


Nonlinrar  acoustic  theory  is  used  to  show  that  the  cutofT  frequencies  and  the  resonant  frequencies  of  I 
modes  in  acoustic  waveguides  of  finite  length  depend  upon  the  mode  amplitude.  The  amplitude  «.  of  I 
mode  n >0,  produced  by  a periodic  piston  motion  of  amplitude  6.,  is  also  determined.  It  is  found  j 
that  at  and  near  a resonant  frequency,  while  at  and  near  the  cutoff  frequency  i 

Away  from  th^  frequencies  e,~8„  as  linear  theory  predicts.  Thus  the  infinite  amplitude,  which  I 

linear  theory  yields  at  cutoff  and  resonance,  is  avoided.  The  results  are  obtained  simply  by  using  the  ' 
result  for  the  propagation  constant  as  a function  of  mode  amplitude  given  by  J.  B.  Keller  and  M.  ' 
H.  Millman  (J.  Acoust.  Soc.  Am.  49,  329-333  (1971)).  i 
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2.  Acoustic  vibrations  in  tubes. 


ARO-28,  "Nonlinear  Resonant  Oscillations  in  Open  Tubes" 
by  Brian  R.  Seymour  and  Michael  P.  Mortell,  J.  Fluid  Mech., 
60,  4,  733-749  (1973). 

A gas  in  a tube,  one  end  of  which  is  open,  is  driven  by  a periodic  applied  velocity 
or  pressure  at  or  near  a resonant  frequency.  Damping  is  introduced  into  the 
system  by  radiation  of  energy  through  the  open  end.  It  is  sho^vn  that  shocks  are 
possible  at  an  open  end  and  that  there  is  a critical  level  of  damping  which 
ensures  a continuous  gas  response  for  all  frequencies.  At  the  critical  level  the 
amplitude  of  the  response  is  0(el),  where  e is  the  amplitude  of  the  input,  and 
it  is  boimded  by  the  amplitude  predicted  by  linear  theory.  There  is  agreement 
with  the  qualitative  experimental  results  available. 


ARO-30,  "Standing  Waves  in  an  Open  Pipe;  A Nonlinear 
Initial-Boundary  Value  Problem"  by  Michael  P.  Mortell, 
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C.  Perturbation  theory. 


1.  Nonlinear  oscillations, 
a.  Nonlinear  wave  equation. 

ARO-12,  "Asymptotic  Solutions  of  Initial  Value  Problems 
for  Nonlinear  Partial  Differential  Equations"  by  Joseph  B. 

Keller  and  Stanley  Kogelman,  SIAll  J.  Appl.  Math.,  4,  748-758 

(1970)  . 


ARO-39,  "Periodic  Solutions  of  Nonlinear  Wave  Equations  in 
N-Dimensional  Space"  by  L.  Ting,  SIAM  J.  Appl,  Math.,  in  press. 

The  solvability  conditions  for  a linear  inhomogeneous  wave 
equation  with  solutions  periodic  in  space  and  time  are  established 
as  integral  conditions.  These  conditions  are  incorporated  in  a 
systematic  perturbation  theory  for  the  determination  of  the 
periodic  solutions  of  the  nonlinear  wave  equation  u^^  - Au 
= ef(u)  with  e as  the  small  parameter.  Applications  to  boundary 
value  problems  are  presented. 


ARO-44,  "Slowly  Modulated  Oscillations  in  Nonlinear  Diffusion 
Processes"  by  Donald  S.  Cohen,  Frank  C.  Hoppensteadt  and 
Robert  M.  Miura,  SIAM  J.  Appl.  Math.,  33,  2,  217-229  (1977). 


It  IS  shown  here  that  certain  systems  of  nonlinear  (parabolic)  reaction-diffusion 
equations  have  solutions  which  are  approximated  by  oscillatory  functions  in  the  form  /?(f 
s^re  P{t  ) represents  a sinusoidal  oscillation  on  a fast  time  scale  t*  and  R((-ct)  represents  a 
stowly^arymg  modulating  amplitude  on  slow  space  (f)  and  slow  time  (r)  scales.  Such  solutions 
desenbe  phcTOmena  in  chemical  reactors,  chemical  and  biological  reactions,  and  in  other  media  where 
a suble  osallauon  at  each  point  (or  site)  undergoes  a slow  amplitude  change  due  to  diffusion 


mf 
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b.  Nonlinear  ordinary  differential  equations. 

I 

ARO-1,  "Periodic  Vibrations  of  Systems  Governed  by  Non- 
linear Partial  Differential  Equations"  by  Joseph  B.  Keller 
and  Lu  Ting,  Conun.  Pure  Appl.  Math.,  1^,  4,  371-420  (1966). 

ARO-41,  "Iterated  Averaging  Methods  for  Systems  of  Ordinary 
Differential  Equations  with  a Small  Parameter"  by  S.  Persek  and 
F.  C.  Hoppensteadt , Comm.  Pure  Appl.  Math.,  in  press.  Exten- 
sions of  Bogoliuboff ' s averaging  methods  are  derived  by  multi-  | 

time  perturbation  procedures.  These  procedures  are  shown  to 
provide  accurate  approximate  solutions  in  cases  where  the 

averaged  system  has  stationary  and  decaying  modes.  Estimates  of 

2 

the  approximations  are  found  on  large  0{l/e  ) and  infinite  time 
intervals,  depending  on  stability  properties  of  certain  second 
order  averaged  systems.  Problems  not  amenable  to  Bogoliuboff 's 
methods,  such  as  the  construction  of  periodic  solutions  and 
their  transients  near  Hopf  bifurcations  in  quadratic  systems, 
can  be  solved  by  this  iterated  averaging  method. 
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c.  Thermal  convection. 

I 

ARO-2,  "Periodic  Oscillations  in  a Model  of  Thermal 
Convection"  by  Joseph  B.  Keller,  J.  Fluid  Mech.,  3,  599-606 

(1966)  . 

Periodic  oscilhil  ion^  c.'.c  I'oiukI  iii  ti  oao-diincu-kit'a!  model  of  f hormal  eom-cclion. 

I'lio  model  coii.'Is!.->  of  ;i  <l(iid-filic(l  tiilic  bool  into  reeiun^^idar .‘•Itnpo  ntid  .-:l.i(idim.'  I 

in  a vertical  iilaiic.  'J'ho  lliiid  is  heated  at  the  eentre  of  Iho  lower  lioii/mit.d 
segmirnl  and  cooled  at  the  centre  of  the.  ii])j)er  liori/.’.';it;!l  segmenl.  'When  a 

certain  parameter  exceeds  unity,  a periodic  motion  of  the  linid  is  found  in  wldoh  j I 

the  flow  is  always  in  llie  same  I'.iicction  but  in  wliie!'.  die  .‘niecd  vatic--.  In.-riia 
is  iitiimoortant  for  tliis  oseillatioti,  whi,'h  drpend.s  njion  ihe  interplay  hetwe'-ii 
frictional  ai.d  bncAai'icy  fori-es. 
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Nonlinear  difference  and  integral  equations. 


ARO-40,  "A  Nonlinear  Renewal  Equation  with  Periodic  and 
Chaotic  Solutions"  by  Frank  C.  Hoppensteadt , SIAM-AMS  Proceedings 
51-60  (1976). 

A nonlinear  renewal  equation  which  arises  in  several  areas  of  nathe- 
matical  population  theory  is  studied  by  a ccnbination  of  rrathcriatical  and 
numerical  analysis.  The  model  is  characterised  by  two  parameters;  m,  a 
measure  of  the  population's  viability  and  fertility,  and  t,  the  (normalized! 
length  of  the  population's  reproductive  window.  Solutions  .ire  described  for 
all  values  of  these  parameters,  0 £ m < 4 , C £ t < 1 . by  a corliination  of 
multi-time  perturbation  analysis  and  numerical  solution.  In  various  te-  j 

gions,  the  solutions  are  shown  to  be  described  by  Burgers'  fcC,uation,  a 
Korteweg-deVries  equation  and  by  a nonlinear  difference  equation.  Kuineri.-al 
methods  are  used  to  investigate  the  remaining  regions.  1 
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2 . Nonlinear  wave  propagation. 

ARO-10,  "Perturbation  Theory  of  Nonlinear  Electromagnetic 
Wave  Propagation"  by  Joseph  B.  Keller  and  Martin  H.  Millman, 
The  Physical  Review,  181 , 5,  1730-1747  (1969). 


A perturbation  method,  previously  employed  to  treat  nonlinear  boundary-value  problems  involving 
partial  differential  equations,  is  used  to  study  nonlinear  electromagnetic  wave  propagation.  The  problems 
considered  are  I.  electromagnetic  wave  propagation  in  a waveguide  containing  a nonlinear  isotropic  me- 
dium; II.  electromagnetic  oscillations  in  a cavity  containing  a nonlinear  dispersive  anisotropic  medium; 

III.  propagation  and  interaction  of  electromagnetic  w.'ves  in  a nonlinear  dispersive  anisotropic  medium; 

IV.  reflection  of  an  electromagnetic  wave  from  a nonlinear  dispersive  uniaxial  medium.  The  method  avoids 
secular  terms,  which  arise  in  some  perturbation  treatments.  The  results  show  how  the  wave  number,  projta- 
gation  velocity,  and  angle  of  refraction  of  a wave  depend  upon  its  amplitude. 


ARO-24,  "Asymptotic  Theory  of  Nonlinear  Wave 
by  Stanley  Kogelman  and  Joseph  B.  Keller,  SIAM  J. 


Propagation" 
Appl.  Math. , 


24,  3,  352-361  (1973). 


By  extending  the  asymptotic  method  of  G.  E.  Kuzmak.  J.  C.  Luke  has  shown  how  G.  B. 
Whitham's  theory  of  nonlinear  wave  propagation  can  be  derived  directly  from  the  partial  differential 
equation  without  using  the  variational  principle,  in  special  cases.  We  apply  the  same  method  to  a 
more  general  class  of  nonlinear  second  order  partial  differential  equations  or  systems  of  first  order 
equations  containing  a small  parameter  c and  obtain  asymptotic  expansions  of  the  solutions. 


ARO-3,  "A  Modified  Stationary  Principle  for  Nonlinear  Waves" 
by  Frederic  Bisshopp,  J.  Diff.  Eqns.,  5,  3,  592-605  (1969). 

Two  theorems  related  to  propagation  of  modulated,  nonlinear  waves  which  are 
governed  by  stationarity  of  the  integral  of  a Lagrangian  density  are  derived.  They 
are  then  used  in  the  analysis  of  asymptotic  approximations  which  represent  nearly 
uniform  wave  trains.  A relatively  detailed  treatment  of  a simple  case  is  folloxved  by 
a summary  of  some  results  which  apply  to  the  propagation  of  capillary-gravity  waves 
in  a liquid. 


3.  Nonlinear  stability  theory. 


ARO-6,  "Finite  Amplitude  Effect  on  the  Stability  of  a Jet 
of  Circular  Cross-Section"  by  D.  P.  Wang,  J.  Fluid  Mech.,  2i» 
2,  299-313  (1968) . 


The  efifect  of  finite  amplitude  on  the  stable  and  unstable  states  of  a column  of 
an  ideal  fluid  of  circular  cross-section  under  the  action  of  surface  tension  is 
studied.  The  method  of  solution  is  a formal  extension  of  the  linearized  theory;  j 
it  consists  of  assuming  that  the  exact  solution  may  be  expanded  in  a power  , 
series  of  a small  parameter  characterizing  the  amplitude.  The  calculation  is  ; 
carried  out  to  the  point  whore  the  first  non  trivial  term  of  the  finite  amplitude  ‘ 
effect  is  obtained.  For  the  stable  states,  the  result  shows  that  the  characteristic 
wavelength  of  a disturbance  which  appears  to  be  stationary  with  respect  to  an 
observer  is  decreased  by  the  finite  amplitude  effect.  For  the  unstable  states,  it 
reveals  that  the  growth  rate  depends  not  only  on  the  wavelength  and  the  magni- 
tude but  also  on  the  t\"pe  of  disturbance  imposed  initially.  The  last  result  is  a 
direct  consequence  of  the  fact  that  two  independent  types  of  initial  disturbance, 
the  disturbance  of  the  velocity  field  and  the  disturbance  of  the  free  surface,  may 
be  imposed  simultaneously  on  the  jet. 


ARO-25,  "Multiple  Scaling  Techniques  for  Nonlinear 
Stability  Problems"  by  S.  Kogelman,  Proc.  Conf.  on  Stability 
Theory,  Washington  State  Univ. , 179-194  (1972). 

We  consider  a model  equation  which  preserves  the  basic  features  ' 
of  sane  typical  nonlinear  hydrodynamic  stability  problems.  By  in- 
troducing two  time  scales  we  are  able  to  obtain  asymptotic  solutions 
to  the  initial  value  problem.  Furthermore,  introduction  of  an 
additional  spatial  variable  enables  us  to  account  for  modes  whose 
wave  numbers  are  very  near  the  critical  wave  number  predicted  by 


linear  theory. 


I 


ARO-14,  "Transient  Behavior  of  Unstable  Nonlinear  Systems 


with  Applications  to  the  Benard  and  Taylor  Problems"  by 


Stanley  Kogelman  and  Joseph  B.  Keller,  SIAM  J.  Appl.  Math., 


20,  4,  619-637  (1971)  . 


A nonlinear  initial  value  problem  conliiiniiig  tv\o  parameters,  r.  and  is  considered  for  a 
function  For  each  r,  uU)  is  a vector  in  a Hilbert  space.  When  c = 0 the  problem  has  a steady 
state  solution  for  any  value  of  a.  This  solution  is  assumed  to  be  linearly  stable  for  /.  < and  linearly 
unstable  for  i.  > /.,  for  some  value  This  means  that  the  linear  problem  for  the  derivative  at 
£ = 0 has  e.vponentially  growing  solutions  for  /.  > but  not  for  /.  < The  solution  «(()  is  found  for 
e small  and  /.  slightly  larger  than  It  is  found  that  ii  cimiains  one  unstable  mode  which  initially  in- 
creases exponentially  but  ultimately  approaches  a constant.  v\  hile  all  other  modes  decay,  up  to  terms  of 
order  r.'.  To  this  order  «(r)  appro.iches  a steady  state  solution  id  x ) dilTcre.nt  from  as  / increases. 
The  method  of  analysis  is  the  two-time  perturbation  method.  The  methvsd  and  re.Milts  are  illustrated  by 
applying  them  to  the  Benard  problem  of  conveelive  heat  transfer  through  a horizontal  layer  of  viscous  j 
fluid  heated  fro.m  below,  and  the  Taylor  problem  of  the  instability  of  the  Couette  flow  of  a viscous  fluid 
between  coaxial  rotating  circular  cylinders. 


ARO-27,  "Nonlinear  Stability  Theory"  by  Joseph  B.  Keller, 


Proceed.  Conference  on  Stability  Theory,  Washington  State 


University,  85-98  (1972) . 


ARO-37,  "Nonlinear  Stability  Analysis  of  Static  States 


which  Arise  Through  Bifurcation"  by  F.  Hoppenstead  and  N.  Gordon, 


Comm.  Pure  Appl.  Math.,  28,  355-373  (1975). 


Multi-time  expansion  solutions  are  constructed  for  systems  of  non- 
linear partial  differential  equations  near  bifurcated  states.  The  class 
of  problems  studied  includes  several  problems  of  interest  in  stability 
of  fluid  flow,  such  as  the  Benard  and  Taylor  problems,  problems  of  non- 
linear diffusion,  as  well  as  many  initial  value  problems  for  systems  of 
ordinary  differential  equations.  The  main  result  gives  a rigorous 
derivation  of  an  expansion  for  the  solution  as  a sum  of  expansions  for 
the  bx  .ced  state  and  one  each  for  slow  and  fast  time  transients  to 
that  state.  The  expansion  for  the  transient  gives  an  approximation  to  the 
domains  of  attraction  of  stable  bifurcated  states.  Also,  as  a consequence 
of  the  method,  we  give  a rigorous  derivation  of  the  Landau  equation  for 
these  systems. 


I 
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ARO-29,  "On  the  Asymptotic  Solution  of  a Two-Parameter  I 

Boundary  Value  Problem  of  Chemical  Reactor  Theory"  by  Jye  Chen 
and  R.  E.  O'Malley,  Jr.,  SIAM  J.  Appl,  Math.,  26,  4,  717-729 
(1974)  . 

The  paper  provides  an  asymptotic  solution  to  boundary  value  problems  of  the  fcr~  j 
c,y"  + a(x)y'  = /(x,  y,  z),  t,£,r"  + IKx)z'  = g(x,  y,  z), 
y'(0)  = z'(0)  = 0.  £,y-(l)  + >U)  = £.£,z'(l)  + z(l)  = B, 

where  £|  and  £,  are  small  positive  parameters  simultaneously  tending  to  zero,  a and  b are  strictly 
positive,  and  a,  b.f,  and  g are  infinitely  differentiable  in  their  arguments. 

ARO-34,  "Asymptotic  Stability  in  Singular  Perturbation 
Problems.  II;  Problems  Having  Matched  Asymptotic  Expansion 
Solutions"  by  Frank  Hoppensteadt , J.  Diff.  Eqs.,  ]^,  3, 

510-521  (1974). 

The  stabilit>’  of  systems  of  ordinary  differential  equations  of  the  form  | 

dxldt  --  /(f,  .V,  «).  € dyfdt  ^ g{t,  x,  y.  c), 

%vhere  cis  a real  parameter  near  zero,  is  studied.  It  is  shown  that  if  the  reduced 
problem 

dxldt  = /(r.  X.  y,  0),  0 = g(t,  X,  y,  0), 

is  stable,  and  certain  other  conditions  which  ensure  that  the  method  of  matched 
asymptotic  expansions  can  be  used  to  construct  solutions  are  satisfied,  then  the  | 
full  problem  is  asymptotically  stable  as  ( ->  oc,  and  a domain  of  stability  is 
determined  which  is  independent  of  t.  Moreover,  under  certain  additional 
conditions,  it  is  shown  that  the  solutions  of  the  perturbed  problem  have  limits 
as  ( — * 00.  In  this  case,  it  is  shown  how  these  limits  can  be  calculated  directly 
from  the  equations 

/(«.  X,  y,  c)  = 0 g(oo,  *,  y,  »)  = 0 

as  expansions  in  powers  of  c.  I 
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ARO-5,  "A  Non-Linear  Singular  Perturbation  Problem  Arising 
in  the  Study  of  Chemical  Flow  Reactors"  by  R.  E.  O'Malley,  Jr., 
J.  Inst.  Maths.  Applies.,  6,  12-20  (1969). 


This  paper  presents  a method  of  obtaining  the  complete  asymptotic  solution  of  boundary 
value  problems  of  the  form 

ty"-b(x)y'-g(x,y)  ==-•  0 
ey'(0)  = y(0)-a 

y(i)  = r> 

for  xe  [0, 1]  where  b(x)  is  strictly  positive  and  for  e small  and  positive.  Physically,  the 
problem  arises  in  determining  the  steady-state  concentration  of  a substance  in  a chemical 
flow  reactor.  A “two-variable”  e.xpansion  procedure  is  used. 


ARO-9,  "Perturbation  Theory  of  Nonlinear  Boundary-Value 
Problems"  by  Martin  H.  Millman  and  Joseph  B.  Keller,  J.  Math. 
Phys.,  2,  342-361. 


A systematic  jicrturbation  theory  is  presented  for  the  analysis  of  nonlinear  problems.  The  lowesl-ordcr 
result  is  just  that  obtained  by  linearizing  the  problem,  and  the  higher-order  terms  arc  the  solutions  of 
inhomogenessus  linear  problems.  The  essential  feature  of  the  method  is  the  procedure  for  avoiding  secular 
terms,  which  is  based  on  the  Lindstedt-Poincare  technique  employed  in  celestial  mechanics.  The  method 
is  applied  to  the  following  nonlinear  boundary  value  problems:  (I)  tcm('erature  distribution  due  to  a 
nonlinear  heat  source  or  sink;  (2)  self-sustained  oscillations  of  a system  with  inlinitely  many  degrees  of 
freedom;  (J)  forced  vibrations  of  a "string”  with  a nonlinear  restoring  force;  (4)  superconductivity  in  a 
body  of  arbitrary  shape  with  evternal  magnetic  field;  (5)  superconductivity  in  an  iiirmiie  film  with 
parallel  magnetic  field;  (6)  comparison  of  solutions  of  the  Uartree.  f'oek.  and  Schrbdinger  equations  for 
the  helium  atom.  The  results  in  each  case  are  dilTcrent  both  qualitatively  and  quantitatively  from  those 
of  the  littear  theory. 


ARO-22,  "Justification  of  Matched  Asymptotic  Expansion 
Solutions  for  Some  Singular  Perturbation  Problems"  by 
Frank  Hoppensteadt , Amer.  Math.  Society,  337-341  (1973)  . 

Several  types  of  problcuis  involving  non-linenr  partial  differential  equation.-, 
having  entail  parameters  arc  described  v.’hich  are  amenable  to  the  method  or 
matched  aayniptotlc  expansions.  This  method  gives  approximate  so.Uttion.s  to  Ihtt 
problems  V7hicl\  are.  valid  over  large  time  intervals.  This  theory  is  applic;,b]c 
to  many  problc;ns  involving  the  il.-iv  i er-StoUcs  equal  ions  and  to  quite  gitacr.i] 
operator  cquatjons. 


ARO-4,  "A  Boundary  Problem  for  Certain  Nonlinear  Second 
Order  Differential  Equations  with  a Small  Parameter"  by 
R.  O'Malley,  Arch.  Rat.  Mech.  and  Anal.,  2^,  66-74  (1968). 

ARO-11,  "Asymptotic  Series  Solutions  of  Some  Nonlinear 
Parabolic  Equations  with  a Small  Parameter"  by  F.  Hoppensteadt, 
Arch.  Rat.  Mech.  and  Anal.,  284-298  (1969). 

ARO-17,  Bifurcation  Theory  and  Nonlinear  Eigenvalue 
Problem,  W.  A.  Benjamin  Co.,  New  YorJc  (1969). 

ARO-38,  "Differential  Equations  Having  Rapidly  Changing 
Solutions:  Analytic  Methods  for  Weakly  Nonlinear  Systems" 
by  F.  C.  Hoppensteadt , J.  Diff.  Eqns. , 22,  237-249  (1976). 

Initial  value  problems  for  weakly  nonlinear  systems  of  differential  equations 
are  studied  here.  A perturbation  result  is  obtained  which  is  shown  to  give  a 
valid  approximation  to  the  solution.  The  approximation  is  determined  by 
solving  two  associated  problems,  and  it  is  especially  derived  for  problems 
which  have  rapidly  oscillating  and  rapidly  decaying  solutions. 
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5.  Continuum  mechanics,  etc. 

ARO-7,  "Rossby  Waves  in  the  Presence  of  Random  Currents” 
by  Joseph  B.  Keller  and  George  Veronis,  J.  Geophys.  Res., 

74,  8,  1941-1951  (1969). 


The  effect  of  weak  random  currents  on  time  periodic  Rossby  waves  is  determined  by  a 
perturbation  method  devised  previously.  As  s result  the  dispersion  equation  for  the  aver- 
age Rossby  wave  is  obtained,  and  from  it  the  propagation  velocity  and  decay  rate  are 
found.  These  quantities  are  evaluated  explicitly  for  steady  unidirectional  random  currents. 
From  the  results  simple  forms  are  obtained  when  the  Rossby  wavelength  is  either  large  or 
small  in  comparison  to  the  correlation  length  of  the  current.  For  some  directions  of  propaga- 
tion of  the  wave  the  velocity  is  increased  by  the  current,  whereas  for  other  directions  it  is 
decreased,  and  a physical  explanation  of  this  is  presented.  Random  currents  will  scatter  the 
waves  and  cause  them  to  decay.  If  the  current  motion  can  add  to  the  energy  flux  of  the  waves, 
the  energy  flux  may  grow  drawing  energy  from  the  currents.  It  is  found  that  east-west  currents 
serve  mainly  to  scatter  the  waves  but  north-south  currents  can  contribute  to  the  growth  of 
waves.  The  growth  phenomenon  is  interpreted  as  an  instability  of  the  currents  to  Rossby  waves 
to  which  they  can  give  up  some  of  their  energy.  The  maximum  growth  rate  occurs  for  waves 
traveling  either  northwestward  or  southwestward,  because  the  energy  flux  of  these  waves  is  in 
the  same  direction  as  the  direction  of  the  random  currents. 


ARO-20,  "Flow  of  a Viscous  Fluid  Through  an  Elastic  Tube 
with  Applications  to  Blood  Flow”  by  S.  I.  Rubinow  and 
Joseph  B.  Keller,  J.  Theor.  Biol.,  35,  299-313  (1972). 


It;  st.'ady  flow  of  ;>  \iicoiis  iuiiii  tbiiiuy.h  '.n  clastic  tube  is  dc'.crmfr.,'  1 
tlv.'ci-.'.ically.  It  is  u'.  .uiT.cii  th.ai  I’oi^.uil'c's  liisv  I'OlJv  l<v.a!ly  and  ihal  i.tl 
any  location  tl'.c  r.a.lii>  of  the  lube  is  dcicnviincd  by  llic  tr.'.r.S!'.'c..:l  i 
pr'.js'.irc  ilitTcicncc.  T r.c  stros;,  strain  cm  vo  of  the  tube  iiiaicri.d  is  ; i-in',-.  -J  ‘ 
to  be  known.  Ti'.c  iVix  ihrcay.h  the  tohe  is  found  .ns  a function  of  inict 
ptc-u  o,  outlet  piVvr,'.i!C.  cNtcrf.al  pressure,  tube  Icnylh,  fluid  si^cc^itv 
pnd  ilv’  clastic  p.'op.i  li  j<  of  the  tube.  'I  he  results  arc  di'oloycd  graphic  iMs 
for  the  hi’ni.'m  c.\ter!i.-il  ilia...  artery,  'ilicy  arc  in  jtood  agreement  v,  iiu 
c\p.  ri.'i.eni.al  ste.t.ly  r..ns  I’lortsiireincnts  on  rubber  tubes,  on  veno'.is  reti  i n 
to  the  heart,  and  on  isolated  c.tts’  Imigs.  The  ciTect  of  tube  coll.tpse  is  also 
considered. 
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ARO-21,  "Post  Buckling  Behavior  of  Elastic  Tubes  and  Rings 
with  Opposite  Sides  in  Contact"  by  Joseph  E.  Flaherty, 


Joseph  B.  Keller  and  S.  I.  Rubinow,  SIAM  J.  Appl . Math.,  23,  4, 


446-455  (1972). 


An  elastic  tube  of  circular  cross  section  or  an  clastic  ring  can  buckle  if  p,  the  outside 
pressure  minus  the  inside  pressure,  exceeds  the  buekling  pressure  p^,  ■ As  p increases  above  p^,  the 
cross  section  becomes  somewhat  elliptical.  Ultimately  opposite  sides  touch  at  one  point  at  the  contact 
pressure  . .As  p increases  above  p,,  the  curvature  of  the  cross  section  at  the  points  of  contact  decreases 
until  it  becomes  zero  at  p„..  For  p > p„,  contact  occurs  along  a straight-line  segment  which  increases 
ill  length  as  p increases  The  pressures  p^, . p..,  and  p„,  are  determined  numerically,  and  the  shape  of  the 
cross  section  is  found  for  various  values  of  p Graphs  of  the  results  are  show  n.  Above  Pu.  a similarity 
solution  is  found.  Analogous  results  are  obtained  for  the  nth  buekling  mode,  n > 2.  The  flux  of  an 
incompressible  v iscous  fluid  through  the  buckled  tube  is  also  determined  us  a function  of  p.  The  results 
may  be  applicable  to  the  collapse  of  veins  and  the  flowi  of  blood  through  them. 


I 
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ARO-23,  "Contact’  Problems  Involving  a Buckled  Elastica" 
by  Joseph  E.  Flaherty  and  Joseph  B.  Keller,  SIAM  J.  Appl.  Math., 
2,  215-225  (1973)  . 

The  classical  problem  of  the  buckling  of  a slender  rod  Ian  elastica)  is  reexaminetl 
numerically  and  analytically  taking  into  account  contact  between  different  parts  of  the  rsxl.  The  lowest 
mode  of  the  clamped  elastica  and  the  third  mode  of  the  pinned  elastica  are  treated.  Other  mixles  can 
be  treated  in  a similar  way.  It  is  found  that  for  each  of  these  moiles  there  is  a compressive  buckling 
load  P,  > 0 and  a contact  load  P,  > P,.  The  rod  remains  straight  when  P S P».  the  rod  buckles 
without  contact  for  P,  < P < P,.  and  two  points  of  the  rod  are  in  contact  for  P 2 P,.  The  contact  | 
points  move  tow  ard  the  ends  of  the  rod  as  P increases.  In  each  case  an  asymptotic  formula  is  obtained  i 
for  the  location  of  the  contact  points  as  P - x . In  both  cases  a similarity  solution  is  found  for  the 
shape  of  the  elastica  between  the  two  contact  points.  In  the  pinned  case  there  is  a range  of  loads  just 
above  P,  for  w hich  the  distance  between  the  ends  of  the  rod  increases  us  the  load  increases.  This  indicates 
instability  and  the  possibility  of  “snap  through  ” | 


ARO-31,  "Resonant  Surface  Waves"  by  J.  R.  Ockendon  and 
H.  Ockendon,  J.  Fluid  Mech.,  2,  397-413  (1973). 

Stnall  atiiivlii  lulc  ftu  tvil  liovizoiital  tir  voi  tionl  o.sfillat  ions  of  a contaiiiof  of  liijuid 
ytilli  a f:’rc.sii!f;>oc  tiny  oivo  li.sp  (o  juolioji.s  in  the  )i<|iii<I  on  a .seal;-  tatK-lt  j:r<  itoj 
titan  till-  fo!ci:;o  ninplil tido.  .mi.Ii  sitiialions  arc  analysnt  ami.  in  lim-o 

ca.viv'.^  vhert'  (lie  ifsjion'i.’  is  slill  .siii.i!l  foinjiaro.l  willi  flu;  tliiniaiMdiis  i,f  il,i. 
contiiinof.  c.xiilicil  a-'Viiiiitui in  .suluf  imts  for  the  liijitiil  inotiini  ate  uliiaineil. 


ARO-19,  "Temperature  of  a Nonlinearly  Radiating  Semi- 


Infinite  Solid"  by  Joseph  B.  Keller  and  W.  E.  Olmstead,  Q.  App. 
Math.,  559-566  (1972). 


ARO-43,  'Applications  of  Multi-Time  Methods  to  Pattern 

i 

Recognition  and  Other  Problems"  by  F.  C.  Hoppensteadt  and 
W..  L.  Miranker,  Proc.  Journees  Math,  sur  les  Perturbation 
Singulieres  (1976) , Springer-Verlag,  in  press. 

Systems  of  difference  equations  containing  small  parameters  are 
studied  by  a constructive  perturbation  scheme  analogous  to  the  one 
developed  by  the  authors  for  the  study  of  differential  equations. 

The  method  results  in  an  averaging  procedure  for  difference  equations, 
and  it  is  particularly  well  suited  to  certain  highly  oscillatory,  non- 
linear systems.  The  method  is  applied  to  problems  for  population 
genetics,  pattern  recognition,  regression  analysis,  and  the  numerical 
analysis  of  still  differential  equations. 


ARO-26,  "Some  Trends  in  Applied  Mechanics"  by  J.  B.  Keller, 


Proc.  of  the  1st  Iranian  Congress  of  Civil  Eng.  and  Mech., 
May  5-9,  103-124  (1972). 


ARO-15,  "Pulse  Propagation  in  a Laterally  Heterogeneous 
Solid  Elastic  Sphere"  by  Z.  S.  Alterman,  J.  Aboudi  and 
F.  C.  Karal,  Geophys.  J.  R.  Astro.  Soc.,  243-260  (1970). 

The  motion  of  a solid  clastic  sphere  caused  by  an  impuhiec  point  source 
has  been  calculated  by  a finite  diricrei’ce  sclvenic.  Results  arc  obtained 
for  a radially  and  laterally  hctcror'cncous  sphere.  Both  discontinuous 
changes  in  clastic  parameters  and  density,  and  their  continuous  variation 
inside  the  sphere  are  treated.  The  results  sb.ou'  reneclion  and  difiraciion 
effects.  A st.abiliiy  analysis  for  tlic  finite  difference  scheme  in  spherical 
co-ordinates  is  given.  Tltc  iv.eth.od  of  soluiion  presented  here  is  r.cneral 
enough  to  include  ih.c  actual  variations  in  density  and  elastic  parameters 
in  the  liarth.  This  includes  the  \\cll-!;no\vn  models  of  ladial  distribution 
and  the  recently  observed  lateral  inhomogeneity  in  the  mantle. 


ARO-42,  "Quantitative  Theories  of  Carcinogenesis"  by 
Alice  Whittemore  and  Joseph  B.  Keller,  SIAM  Review,  in  press. 

Various  mathematical  theories  of  the  genesis  of  cancer  are 
presented.  They  involve  the  transformation  of  a normal  cell  and  the 
subsequent  proliferation  of  the  transformed  cell  and  its  descendants 
to  form  a tumor.  The  theories  predict  the  rate  of  tumor  occurrence 
as  a function  of  time  and  of  the  concentration  of  carcinogenic 
agents,  which  cause  transformation  and  may  also  accelerate  growth. 

The  consequences  of  the  theories  are  deduced  and  compared  with  data 
on  human  cancer  incidence  and  with  data  frc»n  experiments  on  laboratory 
animals.  A new  theory  is  presented  which  includes  features  of  most 
previous  theories,  but  which  goes  beyong  them  in  being  capable  of 
accounting  for  the  results  of  initiation-promotion  experiments. 


ARO-46,  "Peeling,  Slipping  and  Cracking — Some  One-Dimensional 
Free  Boundary  Problems  in  Mechanics"  by  R.  Burridge  and  J.  B.  Keller, 
SIAM  Review,  in  press. 

We  examine  several  problems  involving  a one-dimensional  continuum  with 
a free  boundary.  Apart  from  their  intrinsic  interest,  they  illuminate  the 
related  but  far  more  complex  problems  which  arise  in  fracture  mechanics,  in 
frictional  sliding  such  as  may  occur  in  shallow  earthquakes,  and  in  the 
testing  of  adhesives. 

First  we  treat  peeling  of  adhesive  tape  from  a surface.  The  tape  is 
regarded  as  a flexible  string,  part  of  which  adheres  to  the  surface  and 
part  of  which  is  separated  from  it.  The  free  boundary  is  the  point  (or 
points)  at  which  separation  takes  place  and  the  problem  is  to  find  the 
trajectory  of  this  point  (or  points)  together  with  the  motion  of  the 
separated  section  of  the  string.  We  treat  both  infinitesimal  and  finite 
deformations. 

! 
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We  also  treat  the  related  problem  of  the  sliding  of  a segment  of  a 


stretched  string  pressed  against  a rough  plane  under  the  influence  of 
transverse  loads.  Adhesive  forces  may  be  either  present  or  absent  at 
the  end  points  of  the  moving  segment,  which  form  the  free  boundaries. 
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II . Summary. 

The  research  projects  supported  by  these  grants  include 
investigations  into  problems  of  nolinear  wave  propagation, 
nonlinear  electrical  circuits,  propagation  in  wave  guides, 
thermal  convection  and  continuum  mechanics.  The  success  of 
these  projects  rested  on  our  unified  and  coherent  development 
and  applications  of  multi-scale  perturbation  methods  to  study 
the  complicated  nonlinear  systems  which  describe  these  physical 
phenomena. 
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